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Abstract. Determining the conditions under which a black hole can be produced is a long- 
standing and fundamental problem in general relativity. We use numerical simulations of 
colliding self gravitating fluid objects to study the conditions of black-hole formation when the 
objects are boosted to ultrarelativistic speeds. Expanding on previous work, we show that the 
collision is characterized by a type-I critical behaviour, with a black hole being produced for 
masses above a critical value, M c , and a partially bound object for masses below the critical 
one. More importantly, we show for the first time that the critical mass varies with the initial 
effective Lorentz factor (7) following a simple scaling of the type M c ~ if (7)" 1,0 , thus 
indicating that a black hole of infinitesimal mass is produced in the limit of a diverging Lorentz 
factor. Furthermore, because a scaling is present also in terms of the initial stellar compactness, 
we provide a condition for black-hole formation in the spirit of the hoop conjecture. 

PACS numbers: 04.25.Dm, 04.25.dk, 04.40.Dg, 95.30.Sf, 97.60.Jd97.60.Lf 
1. Introduction 

There is little doubt that black-hole formation represents one of the most intriguing and 
fascinating predictions of classical general relativity. There is abundant astronomical evidence 
that black holes exist, and a number of considerations supporting the idea that under 
suitable conditions gravitational collapse is inevitable Q. In addition, there is overwhelming 
numerical evidence that black-hole formation does take place in a variety of environments EB. 
Yet, a rigorous definition of the sufficient conditions for black-hole formation is still lacking. 
Hence, it is not possible to predict whether the collision of two compact objects, either stars 
or elementary particles, will lead to the formation of a black hole. 

However, we are not without ideas and the intuitive hoop conjecture proposed by Thorne 
in the '70s, provides us with some reasonable guidelines [3 J. We recall that the conjecture 
states that a black hole is formed if an amount of "mass-energy" E can be compressed to 
fit within a hoop with radius equal or smaller than the corresponding Schwarzschild radius, 
i.e. if i?hoo P < R s = 2GE/c 4 , where G is gravitational constant and c the speed of light. Even 
though it can be made precise under particular circumstances [4], the hoop conjecture is not 
meant to be a precise mathematical statement and, in fact, it is difficult to predict if the above- 
mentioned collision will compress matter sufficiently to fit within the limiting hoop. Loosely 
speaking, what is difficult is to determine which part of the "kinetic energy" of the system can 
be accounted to fit within the hoop. Since at the collision the conversion of kinetic energy 
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into internal energy is a highly nonlinear process, any quantitative prediction becomes rapidly 
inaccurate as the speeds involved approach that of light. 

As stated above, the hoop conjecture is purely classical. A quantum-mechanical 
equivalent is not difficult to formulate, although not very stringent, as it simply implies that 
a black hole will be formed at Planck-energy scales. The predicting power does not improve 
significantly when considering the conditions of black-hole formation in higher-dimensional 
theories of gravity (see, e.g. J5l|6l[7]]). In these frameworks, the energy required for black-hole 
formation might be significantly smaller [5], thus providing the possibility of producing them 
in the Large Hadron Collider (LHC) (8|, but no firm conclusion has been reached yet. 

Clearly, although numerical simulations represent a realistic route to shed some light on 
this issue (see, e.g. , J9l [TOl [TTJ), even the simplest scenario of the collision of two compact 
objects at ultrarelativistic speeds is far from being simple and it is actually very challenging. 
A first step was taken by Eardley and Giddings [12], who have studied the formation of a 
black hole from the head-on collision of two plane-fronted gravitational waves with nonzero 
impact parameter (previous work of D'Eath and Payne [131 EH UM using different methods 
had considered a zero impact parameter). In all of these analyses each incoming particle 
is modelled as a point particle accompanied by a plane-fronted gravitational shock wave 
corresponding to the Lorentz-contracted longitudinal gravitational field of the particle. At 
the instant of collision the two shock waves pass through one another and interact through a 
nonlinear focusing and shearing. As a result of their investigation, a lower bound was set on 
the cross-section for black hole production, i.e. a > 32. 5 (GE /2c 4 ) 2 , where E is the centre- 
of-mass (lab) energy. More recently, and in a framework which is closer to the one considered 
here, this problem has been investigated by Choptuik and Pretorius lfT6l . who studied the 
collision of two classical spherical solitons, with a total energy of the system in the lab frame 
E = 27 b m c 2 , where m is the "rest-mass", 7 b = 1/ y/l — v^/c 2 and v h the boost velocity. 
They were then able to show that for collisions with sufficiently high boost, i.e. j h > 2.9, a 
black hole can be formed. 

In this work we report for the first time on black-hole production from the collision of two 
compact, self-gravitating, fluid objects boosted at ultrarelativistic speed^f} There are several 
important differences with the previous investigations in [131 [HI Q31 Q21 [16]]. Differently 
from [131 [EH [EH [13, in fact, our colliding objects are not in vacuum and are not treated 
as point particles. Rather, our relativistic stars are obviously extended and self-gravitating 
objects, thus with a behaviour that is intrinsically different. Also, differently from lfT6ll . our 
objects are not described as scalar fields, but as fluids and thus represent a more realistic 
description of baryonic matter, such as the one employed when simulating relativistic heavy- 
ion collisions |[T8ll . These intrinsic differences also make the comparison with the works 
of H3l [HI [151 [13 very hard if possible at all. On the other hand, many analogies exist with 
the collision of bosons stars considered in lfT6l . and that, as we will discuss below, can be 
interpreted within the more general description of black-hole production from ultrarelativistic 
collisions that we provide in this work. 

I While this work was being reviewed, a similar investigation by East and Pretorius has also been reported ifTTl . 
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Table 1. Initial properties of the isolated stars leading to supercritical models nearer to the 
critical line. Reported in the various columns are the central density p c , the gravitational mass 
M, the rest mass M , the radius R, the compactness M/R, the initial velocity and Lorentz 
factor and 75, and the corresponding effective value (7) in units where c = 1 = M . 
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Given these basic differences, it is not surprising that our results also show qualitative 
differences with respect to previous investigations. The most important of these differences 
is that we find that a black hole can be produced even from zero initial velocities if the initial 
masses are large enough; this behaviour is clearly absent in all previous results, where instead 
a critical initial boost is necessary [[131 031 Q21 Q21 O. In addition, we show that for each value 
of the effective Lorentz factor, (7), a critical initial mass exists, M c , above which a black hole 
is formed and below which matter, at least in part, self gravitates. More importantly, both M c 
follows a simple scaling with (7). 

2. Methodology 

The numerical setup employed in our simulations is the same presented extensively in our 
previous works [[191 EOl EEL and the interested reader can find in these references all the 
needed technical details. It is sufficient to remind here that we use an axisymmetric code 
to solve in two spatial dimensions, (x,z), the set of the Einstein and of the relativistic- 
hydrodynamic equations |[T9l . The axisymmetry of the spacetime is imposed exploiting the 
"cartoon" technique, while the hydrodynamics equations are written explicitly in cylindrical 
coordinates. All the simulations use an ideal-gas EOS, p = (T — l)pe, where p is the rest- 
mass density, e the specific internal energy, and T = 2. The initial configurations consist 
of spherical stars, constructed as in [|20l ED after specifying the central density, p c , where 
the latter also serves as parameter to determine the critical model. The stars have an initial 
separation D and are boosted along the z-direction via a Lorentz transformation with boost 
v h /c. To limit the initial violation in the constraints, D is chosen to be sufficiently large, 
i.e. D = 240 M , and we use an optimal composition of the two isolated- star solutions that 
will be presented in a longer paper. The grid has uniform spacing A = 0.08(0.06) M with 
extents x/M Q e [0,80] andz/M e [0, 150(200)], where the round brackets refer to the more 
demanding high-boost cases. Reflection boundary conditions are applied on the z = plane, 
while outgoing conditions are used elsewhere. A summary of the initial properties of the 
isolated stars leading to supercritical models nearer to the critical line is presented in Table [Tj 
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Figure 1. Representative snapshots of the rest-mass density, p in units where c = 1 = M 
(top row), of the Lorentz factor, 7 (middle row), and of the local fluid energy, — uq (bottom 
row), for subcritical models with an initial small boost v^/c = 0.3 (left panel) or a large one 
Vb/c = 0.8 (right panel). Note that the post-collision flow is essentially jet-like for the low- 
boost case (left panel), while essentially spherical for the high-boost case (right panel); in this 
latter case, most of the matter is unbound. 



3. Basic dynamics 

The dynamics of the process is rather simple. As the two stars approach each other, the initial 
boost velocity increases as a result of the gravitational attraction, leading to a strong shock 
as the surfaces of the stars collide. In the case of supercritical initial data, i.e. of stars with 
masses above a critical value, M c , a black hole is promptly produced and most of the matter is 
accreted. Conversely, in the case of subcritical initial data, i.e. of stars with masses below M c , 
the product of the collision is a hot and extended object with large-amplitude oscillations. Part 
of the stellar matter is unbound and leaves the numerical grid as the product of the collision 
reaches an equilibrium. 

Figure [T] shows snapshots at representative times of the rest-mass density, p (top row), 
of the Lorentz factor, 7 = (1 — v l Vij c 2 ) -1 / 2 (middle row), and of the local fluid energy, 
—uq (bottom row), for two subcritical models. The left panel, in particular, refers to a 
binary boosted at v h /c = 0.3. Note that the stars are strongly compressed by the collision, 
with the rest-mass density increasing exponentially. The merged object expand in a jet-like 
fashion along the z-direction, with the bulk of the matter being accelerated up to 7 ~ 16, or 
equivalently, v/c ~ 0.998, but then settling on much slower flows with 7 < 2.1. Furthermore, 
the front of the jet has — u > 1 indicating that part of the shocked matter has sufficient energy 
to have become gravitationally unbound. As a result, the rest-mass density at the center of 
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the merged object is smaller than the maximum density of the initial configuration, although 
the origin still represents the region where the density is the largest. The right panel, on 
the other hand, refers to a highly-boosted binary, i.e. with v h /c = 0.8, with each star being 
initially highly distorted by the Lorentz contraction. Also in this case, the stars are strongly 
compressed by the collision, but the merged object expands in a spherical blast- wave fashion, 
with an almost spherical distribution of matter and bulk 7-factor. The latter reaches values as 
large as 7 ^ 30, or equivalently, v/c ~ 0.999, which, in contrast with the low-boost case, do 
not decrease in time. It is worth emphasizing that these are the first calculations exploring such 
ultrarelativistic regimes in strong-curvature fields. As a comparison, the typical bulk Lorentz 
factors obtained in the merger of binary neutron stars in quasi circular orbits is 7 ~ 1.03 [22]. 
The very large kinetic energies involved in the collision are sufficient to make a very large 
portion of the stellar matter unbound, as clearly shown by the bottom-right panel of Fig. [TJ 
which reports the local fluid energy. The rest-mass density distribution in the expanding 
blast wave has a minimum at the origin, where a large rarefaction is produced by the matter 
expanding as an ultrarelativistic thick shell. 

The marked transition from a jet-like outflow, not too dissimilar from the simple Bjorken 
flow used to model the very early states of relativistic ion-collisions E3ll . to a shell-like 
structure, not too dissimilar from "transverse expansion" modelled in the subsequent stages of 
relativistic ion-collisions (see G4ll and references therein), signals that it is not unreasonable 
to extrapolate some of the results presented here also to the collision of ultrarelativistic 
elementary particles. 

The transition from the two qualitatively-different regimes discussed above is further 
confirmed by the evolution of the rest-mass normalized to the initial value M , which is shown 
in Fig. [2] for the low (top panel) and high-boost (bottom panel) cases. More specifically, we 
show the fraction of the matter for which the local fluid energy — uq > 1 and which we 
therefore consider to be unbound, M ub /M (blue dashed line). We recall that —uq > 1 is 
a necessary but not sufficient condition for matter to be unbound, valid only in stationary 
spacetimes and ours never really are. However, this condition is probably also the only way 
to measure fraction of unbound matter. Also shown is the corresponding bound fraction 
M h /M (red solid line), defined as the complement to the total rest-mass M (black solid 
line), which is conserved essentially to machine precision. Note that this measurement is 
exact only in an axisymmetric and stationary spacetime, and therefore the estimates of M ub 
and M ub are sensible only in their asymptotic values. Clearly, the unbound fraction is just a 
few percent of the total rest-mass in the case of a low-boost collision, with most of the matter 
being confined in the selfgravitating "star". This is to be contrasted with what happens for 
a high-boost collision, where the unbound fraction is ~ 100% of the total rest-mass. This 
behaviour provides a strong indication that, at least for subcritical collisions, the role played 
by gravitational forces is a minor one as the kinetic energy is increased. This is what happens 
in the collision of two particles at ultrarelativistic speeds, where all of the matter is obviously 
unbound. 
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Figure 2. Evolution in collision-retarded time of the fraction of the matter which is unbound, 
M u b/M (blue dashed line), and the corresponding bound fraction M^/Mq (red solid line), 
defined as the complement to the total rest-mass M (black solid line). 

4. Critical behaviour and scaling 

A remarkable property of the head-on collision of compact stars is the existence of type-I 
critical behaviour, which was first pointed out in ll25ll and subsequently extended in [|T9l . In 
essence, in these works it was found that when considering stars with initial zero velocity 
at infinity, it is possible to fine-tune the initial central density p c (and hence the mass) near 
a critical value, p*, so that stars with p c > p* would collapse eventually to a black hole, 
while the models with p c < p* would eventually lead to a stable stellar configuration. As 
a result, the head-on collision of two neutron stars near the critical threshold can be seen 
as a transition in the space of configurations from an initial stable solution over to a critical 
metastable one which can either migrate to a stable solution or collapse to a black hole ETTl . 
As the critical limit is approached, the survival time of the metastable object, r eq , increases as 
r eq = -A In \p c - p% with A - 10 El ED. 

Although the free-fall velocities considered in [[25l [TH were very small, we find that 
the critical behaviour continues to hold also when the stars are boosted to ultrarelativistic 
velocities and that the threshold p* is a function of the initial effective boost. It should be 
emphasized that determining p* becomes especially challenging as the 7-factor is increased 
and the dynamics of the matter becomes extremely violent, with very strong shocks and 
rarefaction waves. However, we were able to determine the threshold for all the range of initial 
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boosts considered, i.e. v^/c e [0, 0.95], % e [1, 3.2], and even to a reasonable accuracy, 
e.g. pi = (3.288023 ± 0.000003) x 10 14 g/cm 3 , for the initial boost of v h /c = 0.3. 

The existence of critical behaviour near which the details of the initial conditions become 
irrelevant and which is the same at different boosts, i.e. A does not depend on 7 nor on 
Pc ([I2H [26]] have shown there is "universality" when varying 7 and fixing p c ), gives us 
a wonderful tool to explore the conditions of black-hole formation also far away from the 
masses and 7-factors considered in this paper. This is illustrated in Fig. |3j which reports the 
gravitational mass of the isolated spherical star as a function of the effective initial Lorentz 
factor 

™~ (JdVT„ u nW) > { ) 

where is the stress-energy tensor, n M is the unit normal to the spatial hyperspace with 
proper volume element dV, and the index refers to quantities measured in the initial 
unboosted frame. We note that we introduce the definition of an effective Lorentz factor ([]} 
because the stars are extended and thus the 7-factor will be different in different parts of 
the star. Expression ([T]), on the other hand, can be taken as ratio of the energies measured 
in the boosted and unboosted frames, and hence a generalization of the Lorentz factor for a 
point particle (Indeed (7) —> 1 for v b —> 1). Of course, other parametrizations are possible, 
still leading to scaling laws, but with slightly different exponents. An extended discussion 
of this point will be presented in a longer paper. Filled circles indicate initial data leading 
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to a black hole, while triangles indicate initial data leading to a "star", whereby we mean an 
object which is at least in part self gravitating (orange errorbars provide an approximate upper 
limit of ~ 8% to the error in the measurements). Also indicated as a blue solid line is the 
critical line separating the two regions of black hole and star formation (the latter is shown as 
a shaded region). Clearly, the numerical results provide a tight fit of the critical line with a 
power law 

M c /M = K < 7 >-» « 0.92 <7>- L03 , (2) 

and which represents the most interesting result of this work. 

Expression Q offers itself to a number of considerations. First, it essentially expresses 
the conservation of energy. Second, in the limit of zero initial velocities, (7) —t 1, we obtain 
that M c ~ 0.92 M , so that the corresponding total mass, 2M C , is only ~ 12% larger than 
the maximum mass of the relative spherical-star sequence, i.e. M max = 1.637 M . Third, in 
the opposite limit of (7) — >> 00, expression ^ predicts that the critical mass will go zero. 
This is indeed what one would expect: as the kinetic energy diverges, no room is left for 
self gravitating matter, which will all be ejected but for an infinitesimal amount which will 
go into building the zero-mass critical black hole. Fourthly, ([2]) is also in agreement with the 
results in USEl, whereby one can recognize the black-hole formation as the crossing of the 
critical line when moving to larger 7-factors while keeping the rest-mass constant. Finally, the 
scaling relation ([2]) can be expressed equivalently in terms of the original stellar compactness, 
M/R as 

(M/R) c = K f ( 7 )- n/ « 0.08 (7)" L13 . (3) 

Since Mi ab = (j)M is the mass in the lab frame and R is the largest dimension in that frame 
being the transverse one to the motion, the ratio (M lab /i?) c = K 1 \^) 1 ~ n ' ~ K , ( r y)~ - 13 
provides the condition for the amount of energy that, when confined in a hoop of radius R, 
would lead to a black hole, as in the spirit of the hoop conjecture. 

5. Conclusions 

We have computed the first ultrarelativistic collisions of compact fluid stars leading to black- 
hole formation. The Lorentz factor reached in our relativistic-hydrodynamic simulations in 
strong-field regimes are considerably larger than those ever explored in inspiralling neutron- 
star binaries. We find that the properties of the flow after the collision change with Lorentz 
factor, with most of the matter being ejected in a spherical blast wave for large boosts. As in 
previous work, we find that the collision exhibits a critical behaviour of type I and that this 
persists also as the initial boost is increased. This allows us to use a simple scaling law with 
a transparent physical interpretation and a condition for black-hole formation in the spirit of 
the hoop conjecture. 
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